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Abstract —A code design problem for memory devises with 
restricted state transitions is formulated as a combinatorial 
optimization problem that is called a subgraph domatic partition 
(subDP) problem. If any neighbor set of a given state transition 
graph contains all the colors, then the coloring is said to be 
valid. The goal of a subDP problem is to find a valid coloring 
with the largest number of colors for a subgraph of a given 
directed graph. The number of colors in an optimal valid coloring 
gives the writing capacity of a given state transition graph. The 
subDP problems are computationally hard; it is proved to be 
NP-complete in this paper. One of our main contributions in 
this paper is to show the asymptotic behavior of the writing 
capacity C(G) for sequences of dense bidirectional graphs, that 
is given by C(G) = Q(n/\nn) where n is the number of nodes. 
A probabilistic method called Lovasz local lemma (LLL) plays an 
essential role to derive the asymptotic expression. 

I. Introduction 

Recent advanced memory devises such as flash memory 
and phase change memory (PCM) require appropriate cod¬ 
ing for improving its write efficiency. For example, coding 
schemes for flash memories 0G1 0 0 can improve write 
efficiency so that lifetime of the flash memory is lengthened. 
Several constrained coding schemes suitable for PCM has been 
presented m mo. Rank modulation developed by Jiang et 
al. C3 that encodes a message onto a permutation of n- 
elements produced a number of technically and mathematically 
interesting problems. In the paper by En Gad et al. 0, they 
considered a combinatorial problem to find a maximum size of 
the decomposition of a state space into several dominating sets 
of state transition graphs in order to determine the maximal 
rate of compressed rank modulation. 

A problem to partition nodes of a given undirected graph 
into several non-overlapping dominating sets is called a do¬ 
matic partition problem lfl~3ll lf3l . The domatic partition problem 
has been long known in the field of theoretical computer 
science. The problem can be considered as a special type 
of a node coloring problem for an undirected graph. Solving 
the problem in an exact manner is considered to be com¬ 
putationally intractable since it was proved that the domatic 
partition problem is an NP-hard problem fl4l . Several efficient 
approximation algorithms has been developed and analyzed 0 
for the domatic partition problems. 

Our aim of this paper is to provide theoretical foundation 
to the problems called subgraph domatic partition (SubDP) 
problems from a graph theoretic view point. In fact, the 
motivation of this paper came directly from the paper of En 
Gad et al. |0- In our problem setting, a directed graph G that 


represents allowable state transitions of a memory devise is 
initially given. We also assume that a fixed size message should 
be written in the memory for each state transition. In order to 
construct an appropriate encoding and decoding functions with 
high coding rate, we must solve an SubDP problem on G. A 
SubDP problem is a problem to find both of a subgraph of G 
and its domatic partition with largest number of groups; i.e., 
domatic coloring. Although the SubDP problem is so closely 
related to the standard domatic partition problem, it has not 
been discussed in literatures as far as the authors know. 

Our approach to the SubDP problem is summarized as 
follows. We will firstly give a rigid definition of the writing 
capacity for a directed graph that corresponds to the solution 
of an SubDP problem. This can be considered as an abstraction 
that transforms coding problems to combinatorial problems on 
directed graphs. We then invoke Lovasz local lemma (LLL) to 
derive a non-trivial lower bound on the writing capacity. LLL 
is a well-known probabilistic method m that is often used to 
show an existence of certain configuration satisfying all the 
requirements. The lower bound based on LLL reveals asymp¬ 
totic behaviors of the writing capacities C(G) of bidirectional 
graphs that follows C(G ) = il(n/ Inn) for a sequence of 
dense graphs where n is the number of nodes. Furthermore, 
we will show that C(G) = 0(n) for a sequence of extremely 
dense graphs by using Turdn’s lemma that is a fundamental 
lemma in extremal graph theory. 

II. Preliminaries 

In this section, notation used throughout the paper and 
definition related to the SubDP problems will be given. 

A. Notation 

Let G = (V, E ) be a directed graph. For a node v £ V, 
the outbound degree and inbound degree of v are denoted by 
d ou t{v) and di n (v), respectively. The average outbound degree 
e ou t(G) of G is given by 

tout(G ) = j— 'y ' doutiv'). (1) 

' ' v&V 

Note that e ou t{G) = |.E|/|Vj holds due to the relation 
J2 v <zv d °ut(v) = \E\. The symbols S out and A out represent 
minimum and maximum outbound degree of G, respectively. 
In a similar manner, the symbols and A in are the minimum 
and maximum inbound degree. If a directed graph G = (V. E ) 
satisfies (u,v) £ E (v,u) £ E , then the graph G is said 
to be a bidirectional graph. In this paper, we use the term 


“bidirectional” instead of “undirected” in order to show that 
there are two directed edges between u and v explicitly. 

B. Basic assumptions 

Assume that a directed graph G = (V, E ), called a state 
transition graph, is given. The memory devise D associated to 
the graph G can store any v £ V as its state. The state of D 
can be observed from outside. We can change the state of D 
from s £ V to s' £ V if the directed edge (s, s') is included 

in E. Let Af = {1,2called an message alphabet. 
We assume the following scenario. When a request to write 
a message to £ Ai comes, we are allowed to change the 
state of D only once along with an edge in E. After the state 
transition for writing to, it is required that a written message 
in D must be correctly recovered by reading the state of D. 
An appropriate pair of an encoding and a decoding function 
satisfying these assumptions is required to realize a memory 
system based on the devise D. 

C. Encoding function 

The main task of our encoder is to output a next state s' 
of D from an input pair of a message to and a current state 
s. The important constraint imposed on the encoder is (s, s') 
is in E or s' = s. This means that only a state transition (or 
no transition) consistent with the state transition graph G is 
allowed. The following formal definition of encoding functions 
certainly includes this consistency condition. 

Let G = (V", E) be a non-empty subgraph of G. A function 
<j) : V x M —» V satisfying 

Vs £ V, Vm £ M, ( s , </>(s, to)) £ E or s = <f>(s, to) (2) 

is called an encoding function. It should be remarked that an 
encoding function does not need to use all the states in V. 
An appropriate choice of V C V may increase the writing 
capacity of D defined later. 

An encoding process is described as follows. An encoder 
receives a message to £ M to be stored in D as an input. 
The first task of the encoder is to fetch the current state of D. 
Let s £ V be the current state stored in D. We then apply the 
encoding function to the current state s and the input message 
to and obtain the next state s' = <j>{s, to). The last task of the 
encoder is to change the state of D from s to s'. 

Assume that mi,m 2 ,m 3 ,... £ Af are a message se¬ 
quence that is to be encoded sequentially and that So is the 
initial state of D. In this case, we obtain the state sequence 
so, si = 4 >(sq, toi), S2 = </>(si, TO2), ■ ■ ■ that can be considered 
as a walk trajectory on the state transition graph G. From this 
encoding process and the definition of the encoding function, 
it is evident that the state s, is contained in V for any i(i > 1) 
if the initial state of D is in V. In the following discussion, 
we assume that So is a state in V. 

D. Decoding function 

A current written message in D must be correctly recovered 
when we want to read the stored message. This condition 
naturally leads to the following definition of the decoding 
function corresponding to an encoding function. 


Let f> be an encoding function associated with the subgraph 
G = (V, E) C G. A function ft : V — > M satisfying 

Vs £V,Vm £ M, ft(ft(s,m)) = m (3) 

is called an decoding function. 

A decoding process is simple. A decoder first fetch the state 
of D when a read request comes. Then, the decoder applies 
the decoding function to the retrieved state s and obtain a 
recovered message m = i/>(s). The condition Q guarantees 
that any stored message can be correctly recovered. 

E. Writing capacity 

It is intuitively evident that there is a close relationship be¬ 
tween topology of a given state transition graph and allowable 
size of a message alphabet. In the following definition, we will 
introduce a measure that quantifies how much information can 
be written in a single state transition. 

Definition 1 (Writing capacity): Let G be a directed 
graph. If there exist G C G and a corresponding encoder- 
decoder pair (ft, ft) satisfying both of the encoding condition 
(]2} and the decoding condition ([3} for a given pair (G. Ai), 
then the pair (G,A4) is said to be achievable. The writing 
capacity C(G) is defined by 

C(G) = max{|Ad| : (G, Ad) is achievable}. (4) 

Although the definition of the writing capacity is fairly 
simple, the evaluation of this quantity for a given graph G 
is not trivial in general. It might be natural to consider the 
following problems on the writing capacity; (1) how to evaluate 
the writing capacity and how to design encoder and decoder 
systematically, (2) worst case computational complexity to 
evaluate C(G), (3) general upper and lower bounds, (4) 
finding relationships between the writing capacity and graph 
parameters (e.g., maximum degree, degree profile, etc.). In the 
following sections in this paper, some of these problems are 
to be discussed. 

F. Subgraph domatic partition 

In the following sections, we will give several lower and 
upper bounds on the writing capacity. Before going into 
the discussions on such bounds, it might be beneficial to 
obtain the intuition on the problem by observing the following 
interpretation to a coloring problem. 

Assume that a subgraph G C G is given. An (-coloring 
is a map from V to [1, £], where the notation [a, b] represents 
the set of consecutive integers from a to b. Let c : V —» [1,(] 
be an (-coloring of G. For v £ V, c(v) £ [1,(] represents its 
color. 

The neighbor set of v £ V is defined by 

N(s) = {s}U{s' £V\(s,s')£E}. (5) 

Consistency between a decoding function ft and an encoding 
function ft can be translated into the following condition. If 
an (-coloring c on G satisfies 

VsGC, (J {c(it)} = [l,4 (6) 

vGN(s) 


the coloring is called a valid (.-coloring. 

It is evident that a color put on a node represents a message 
symbol associated to the node via a decoding function; namely, 
ip(v) = c(v). If we have a valid coloring, it is trivial to 
construct an encoder function <p satisfying the consistency 
condition ip(<p(s,m)) = m. On the other hand, a valid triple 
(G, (f>, ip) naturally leads to a valid coloring. 

For given G C G, let 

7 (G) = max{£ | 3 a valid coloring on G with f-colors}. 

(7) 

From the definition of 7 , it should be noted that the writing 
capacity can be recast as 

G(G) = max 7 (G). ( 8 ) 

0#GCG 

Finding the largest ( satisfying the condition ([ 6 } for a fixed 
G is equivalent to computing the domatic number of G: The 
domatic number of a graph G is defined as the maximum 
number of node-disjoint dominating sets in G, and it is obvious 
that each color class in a valid ^-coloring constitues a domi¬ 
nating set of G. Since it is allowed to choose an appropriate 
subgraph G C G in our problem, we call the problem © 
the SubDP problem, which requires to identify the subgraph 
G whose domatic number is the maximum. We also extend 
the terminology of domatic number to its subgraph version. 
The subgraph domatic number is defined as the maximum of 
domatic numbers over all sugraphs. 

Example 1: Let us consider a directed graph G = ( V,E) 
illustrated in Fig[T](Left) that represents the adjacency relation¬ 
ships of vertices of 3-dimensional hypercube. Figure |T| (Right) 
provides an optimal valid coloring (G = G). In this case, we 
can let M = [1,4]. A state trajectory of an encoding process 
can be seen as a walk on vertices along with edges of the 
hypercube. 



Fig. 1. A state transition graph and an optimal valid coloring: A node of 
the left graph corresponds to a vertex of 3-dimensional hypercube. A binary 
3-tuple in a node represents its coordinate. A pair of nodes corresponding to a 
pair of adjacent vertices has bidirectional edge connections. This means that a 
state transition changes the binary state label only by one-bit. (Right) Optimal 
valid coloring for this graph. We can prove that C(G ) = 4 in this case. 

Example 2: Figure [2] shows a non-trivial case. Petersen 
graph is a 3-regular graph with 10-nodes. In this case, we 
can prove that G(G) = 3. 



Fig. 2. An optimal subgraph of Petersen graph is denoted by solid lines. The 
optimal valid coloring with 3-colors is indicated with the numbers form 1 to 

3. 

G. subDP is NP-complete 

It is not difficult to show that the subgraph domatic 
partition problem is computationally hard. 

Theorem 1: The subDP problem is NP-complete even if 
the input graph G is d-regular for any d > 3. 

(Proof) It is trivial that the problem belongs to NP. Thus what 
we have to prove is the NP-hardness of the problem. The proof 
is the reduction from full domatic partition problem for regular 
graphs, which is the problem of deciding whether a given 
regular graph is domatically full, that is, the domatic number of 
a given d-regular graph is d+1, or not. It is known that the full 
domatic partition problem is NP-complete for any d> 3 El. 
The core of the reduction is the that the subgraph domitic 
number of any (/-regular graph G becomes d +1 if and only if 
for any G is domatically full, which is verified by the following 
simple observation: Since any subgraph G C G contains a 
node with degree d — 1 or less, its domatic number cannot be 
greater than d. Thus the sugraph domatic number of G can 
become d + 1 if and only if G = G and the domatic number 
of G is d T 1, which implies that G is domatically full. □ 

It is worth pointing out that the computational complexity 
of the (standard) domatic number belongs to much harder class 
than NP-completeness. The domatic number is not only hard 
to compute the exact value, but also hard to be approximated 
even within a factor better than In n 0 . Since the subgraph do¬ 
matic number appears to include more complex combinatorial 
choices, we conjecture that computing the subgraph domatic 
number also belongs to the class extremely harder than NP- 
completeness. However, we unfortunately failed to prove the 
conjecture, and it is still open. Note that applying the strategy 
used in the proof for the (standard) domatic number is quite 
far from trivial. 

III. Main contributions 

The following two theorems express asymptotic behaviors 
of the writing capacity G(G) for bidirectional graphs. These 
results give us insight on the behavior of G(G) for dense graph 
sequences. 

Theorem 2 (Dense graph sequence): Let 
Gn = (Vni -Era), 71 = 1 , 2, . . . 




be a sequence of bidirectional graphs satisfying the following 
conditions. The number of nodes \V n \ is assumed to be n and 
the number of edges grows as \E\ = cm 2 where a(0 < 7 < 1) 
is a real constant. For such a sequence G n , 

c(G n ) = n (^-) (9) 

holds in the regime where n —>• 00. 


Theorem 3 (Extremely dense graph sequence): Let G n = 
(V n: E n ), n = 1, 2,... be a sequence of bidirectional graphs. 
If the number of edges grows as 


\E n \ =n 2 


1-0 



( 10 ) 


then 

G(G„) = Q(n) 


holds in the regime where n —> 00. 


In the following sections, we will establish required lem¬ 
mas to prove these theorems. 


IV. Bounds on writing capacity 

In this section, upper and lower bounds on C(G ) will be 
presented. 


is bounded from below as 


C(G) > 


S ou t(G) + 1 

1 + In 3 + In 2 + 3 In A ’ 


where A = A out (G) + A in (G). 


( 12 ) 


(Proof) We will put ^-colors (£ £ [1, A out (G) + 1]) to nodes 
of G uniformly randomly. A bad event on v € G is the event 
such that there is a missing color in the neighborhood set N(v). 
Precisely speaking, {c(v)\v £ V} ^ [l,f] holds for a bad event 
on v, where c(v) £ [ 1 , (!) represents the color put on the node 
v. Assume that the probability p(v) is the the probability of 
the bad event. By using the union bound, the probability p(v) 
can be upper bounded as follows: 


p(v) < 


i 

Prob [node with color i ^ 7V(u)] 

i— 1 



< ^exp 


W(v)\ ^ 


(13) 


We used the inequality 1 — x < exp(— x) in the last step in the 
above derivation. From this upper bound on p(v), it is easy to 
derive a uniform upper bound. Since |-/V(u)| > 5 ou t(G ) + 1, 
we have p(v) < p for any v £ G where 


A. Upper bound based on maximum degree 

The following upper bound on the writing capacity is 
simple but is fairly useful to determine the writing capacity 
of several special graph classes. 

Lemma 1 (Degree bound): Suppose that a directed graph 
G = ( V,E ) is given. The writing capacity C(G) satisfies 

O(G) < 1 + A out (G). (11) 

(Proof) Since G = (V, E) is a non-empty subgraph of G, it 
is evident that |jV(i;)| = 1 + d^ ut (v) holds for any v £ V 
and that d^ ut (v) < A ou t(G) for any v £ V. This implies that 
the number of colors of a valid coloring for any V can be 
bounded by 1 + A out ( G ) from above. This fact directly leads 
to the claim of the lemma. □ 


P = (A out (G) + 1) exp + 1 j . (14) 

In order to use LLL, we need to have a upper bound on 
the number of bad events that are not mutually independent. 
Let us denote the number of such events by T(v). The number 
T(v) can be bounded from above in the following way: 

T(v) < d out (v) + d in (v) 

+ ( d ou t(v ) + di n (v))(A out (G) + A i n (G)) 

< 2A 2 = D. 

This is because bad events on v and u are mutually independent 
if the length of the shortest path between the nodes v and u 
in G is greater than 2. 


This lemma clarifies the behavior of C(G ) for a sequence 
of sparse graphs. Let G n = (V n: E n ) be a sequence of a sparse 
graphs that satisfies A out (G n ) = 0(1). The degree bound 
£□1 gives inequality C(G n ) < A out (G n ) + 1 that implies 
C(G) = 0(1) as well. 

B. Lower bound based on Lovdsz local lemma 

In this subsection we will show the following lower bound 
on the writing capacity 0(0). The proof of this lemma is 
based on Lovasz local lemma that is a common technique to 
show the existence of valid configuration that satisfies all the 
requirements. The outline of the following argument follows 
the proof of Lemma 2 of Feige et al. 0 . 

Lemma 2 (LLL lower bound): Let G = (V, E) be a di¬ 
rected graph and G be any subset of G. The writing capacity 


LLL guarantees the existence of a valid coloring if 

ep(D + 1) < 1 (15) 

holds. The lefthand side ep(D + 1) can be evaluated as 

ep(D + 1) = e(A out {G) + l)exp + ^ (D + 1) 

< ex pf- ^ (( ; ) + 1 +l+ln(6A 3 )V 


We will consider the condition that the exponent becomes 
non-positive. From the non-positiveness condition 


5out{G) + l +1 + ln(6A 3) < 0 , 


06) 












we have a sufficient condition: 


p < Sout{G) + 1 

— 1 + ln(6A 3 ) 


(17) 


sequence of induced subgraphs Go 2 Gi 7 G 2 3 
According to this choice, we obtain 

&out (G) > max{A out (G), e out (G)/ 2 } 


2G k . 

( 20 ) 


for existence of a valid coloring. Namely, if l is bounded as 
in the above inequality, the LLL condition ( fl5] i is hold. Since 
the number of colors £ is integer, it is reasonable to let 


due to Lemma [3] 

Combining the condition (ITSb and (l20l) . we immediately 


( 21 ) 


£ = 

8 0 ut(G) + 1 

(18) 0 

max{c5 OKt (G), e out (G)/2} + 1 


1 + ln(6A 3 ) 

V ' l — 

1 -h In 3 + In 2 + 3 ln(A OW i (G)) _ 


because our aim is to maximize l. From LLL, it is evident 
that the condition (fill) is a sufficient condition of existence of 
valid ^-coloring for G. □ 

C. Proof of Theorem [2] 

The proof of Theorem [2] can be divided into two parts; 
in the first part of the proof, we will show that there is an 
induced subgraph of a bidirectional graph G that has minimum 
outbound degree larger than the half of average outbound 
degree of G. Since the minimum outbound degree restricts the 
number of colors in a valid coloring, finding an appropriate 
subgraph of G is crucial to increase C(G). In the second part 
of the proof, random coloring is employed and we shall use 
Lemma[2]to derive a sufficient condition in asymptotic regime. 

We firstly will show that any bidirectional graph has an 
induced subgraph with minimum outbound degree larger than 
half of the average outbound degree. A process to obtain an 
induced subgraph with large minimum outbound degree is 
described in the following definition. 

Definition 2 (Subgraph sequence): Let G = (V, E) be a 
bidirectional graph. A sequence of induced subgraphs of G, 

Gq = G D Gi D G 2 D • • • D G k 

is recursively constructed in the following way. If there exists 
a node v £ Gi satisfying 

^w(Gi) > d%u t (v), (19) 

then Gj+i is obtained by removing both v and the adjacent 
edges from Gi. The terminal subgraph Gk does not contain a 
node satisfying the condition ( fl~9l >. The sequence constructed 
in such a way is said to be appropriate. 

The next lemma states that the terminal subgraph Gk is 
what we want. 

Lemma 3: Let G = (V, E ) be a balanced directed graph. 
Suppose also that an appropriate sequence of induced sub¬ 
graphs of G, Go 3 Gi D G 2 3 ■ ■ ■ D Gk, is given. 
The minimum outbound degree of Gk satisfies the inequality 

5out(G K ) > £out(G)/ 2. 

The proof of this lemma is given in Pfl . 


Note that A = 2A out (G) holds for a bidirectional graph. 

Assume that G n = (V n ,E n ) is a sequence of bidirectional 
graph with \E n \ = an 2 . The lefthand side of (ITU can be lower 
bounded as 

max{e OM t(G n )/2,8 

out (G»)} + 1 

A + 31n(A out ) 

where A = 1 + In 3 + In 2. It is clear that we have C(G) = 
f2(n/lnn) under the condition n —>• 00. This completes the 
proof of Theorem U 

As a byproduct of the LLL lower bound, a polynomial time 
approximation algorithm for evaluating the writing capacity 
with approximation ratio 0(ln A out (G)) can be derived where 
A out (G) is the maximum outbound degree of G. 

D. Lower bounds on writing capacity based on Turdn’s lemma 

A graph with extremely high edge density contains a 
subgraph that have a large writing capacity such as complete 
graphs and bipartite graphs. If a given graph G includes such 
a subgraph, its writing capacity can be lower bounded by the 
writing capacity of the subgraph. The result shown in this 
subsection is based on the Turan’s lemma that is a fundamental 
result in extremal graph theory. 

We firstly start our argument on this subsection with the 
following basic lower bound. 

Lemma 4 (Subgraph lower bound): Let G = {V,E) be a 
directed graph and G be a subgraph of G. We have G(G) > 
G(G). 

(Proof) The subgraph G includes a subgraph G'(G' C G C G) 
with G(G)-color valid coloring. Because the subgraph G' is 
also a subgraph of the graph G, it is clear that G(G) becomes 
a lower bound of G(G). □ 

Turan graph is the undirected graph such that it has the 
largest number of edges among undirected graphs with n- 
vertices which does not contain a complete graph K n as a 
subgraph. The number of edges of Turan graph is denoted by 
f r _i(n). Turan proved the following upper bound on the size 
of the Turan graph. 


( 22 ) 

A + 3 In n 


We are ready to go into the proof of Theorem [2] Let G = 
(V. E) be a bidirectional graph. A subgraph G can be chosen 
arbitrary so that 8 ou t(G ) becomes large. We here take the 
following strategy to choose G: If 8 out (G) > e ou t(G)/ 2, then 
let G = G. Otherwise, namely S ou t{G) < e ou t{G)/2, then let 
G = Gk where Gk is the terminal subgraph in an appropriate 


Lemma 5 (Turan’s lemma): 

t r -i{n) < \n 2 - —(23) 
2 r — 1 

The upper bound of the edge size of Turan graph leads to 
the following lemma that guarantees the existence of a large 














clique for an extremely dense graph. This lemma constitutes (Proof) The writing capacity G(G) can be bounded from above 
a main part of the proof of Theorem [3] as 


Lemma 6: Let G = (V. E) be a bidirectional graph. If 


\E\> 


1 , Icml — 2 

—n — --- 

2 \an\ - 1 


(24) 


holds, then G(G) is lower bounded as G(G) > [cmj where 
ct (0 < a < 1 ) is a real constant. 


The proofs of this lemma and Theorem 0 are omitted. The 
forthcoming extended version of this paper will include these 
proofs and also discussion on a polynomial time approximation 
algorithm for subDP problems. 


We will firstly show that the complete graph K n has the 
writing capacity C{K n ) = n. For each node of K n , we will 
put a distinct color from [l,n]. It is clear that this coloring 
is a valid coloring for K n . Thus, we have C(K n ) > n. 
On the other hand, due to the degree bound, we obtain 
another inequality C(K n ) < n. These two inequalities imply 
C{K n ) = n. 

From Turan’s lemma, if the condition ( l24l > is met, then 
the graph G contains Ky an j as a subgraph. By Lemma [4] we 
finally have G(G) > [an\. □ 


The proof of Theorem [3] is given as follows. The coefficient 
appeared in (l24l > can be bounded from above: 


Lq^J - 2 l i k1 __L 

Lcmj — 1 [an] — 1 an 


(25) 


Due to the assumption on the number of edges, there exists 
a positive number C satisfying \E\ > n 2 /2 (1 — C/n) for 

sufficiently large n. By letting a = 1/C, Lemma [ 6 ] guarantees 
C(G) > an. On the other hand, the degree bound givens 
G(G) = 0(n). Combining these, we immediately obtain 
G(G) = O(n) in the asymptotic regime. 


V. Approximation algorithm 

As proved in Subsection , SubDP problems are NP- 
complete. It is natural to try developing an approximation 
algorithm to solve a subDP problem that is equivalent to the 
evaluation of the writing capacity of a given graph. In section, 
we will show that evaluation of the writing capacity of a bidi¬ 
rectional graph with the approximation ratio 0(ln A out (G)) 
can be done in polynomial time. The key to derive the 
approximation ratio is to use a tight upper bound that can 
be computed in polynomial time with n. Let G be a bidireced 
graph with n nodes. A k-core of G is a connected subgraph G 
with the minimum outbound degree k. It is known that k -core 
can be found in linear time if it exists. The following upper 
bound is based on evaluation of fc-cores of G. 

Lemma 7 (k-core bound): Let G be a bidirectional graph 

and 77 (G) be 77 (G) = max{/c|G has a fc-core}. The writing 
capacity G(G) is upper bounded as 

G(G) < 1 + 77 (G). (26) 


G(G) = 

max 7 (G) 

(27) 


0^GCG 


< 

max [1 + <5 out (G)] < 1 + 77 (G). 

(28) 


0#GCG 


The last inequality 

is due to the definition of 77 (G). 

□ 


The process of our approximation algorithm is summarized 
as follows. At the first stage of the algorithm, we will find a 
subgraph G that satisfies 5 ou t(G) = rj(G). Such a subgraph 
G can be obtained in polynomial time. In the second stage of 
the algorithm, an algorithmic version of LLL with setting 


t 


A 


l + v (G) 

1 + In 3 + In 2 + 3 In A 


(29) 


is exploited. It can find a valid coloring with G-colors in 
polynomial time. Closely related approximation algorithms 
for domatic partition problem is intensively discussed in 0 - 
Combining the upper bound (l26b and the lower bound < [29l) . we 
immediately have the following upper bound on the approx¬ 
imation ratio C(G)/£* = 0(ln A out (G)), which is achieved 
with the above algorithm in polynomial time. 
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